Measurement of the Electric Current in a Kpc-Scale Jet 
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ABSTRACT 

We present radio emission, polarization, and Faraday rotation maps of the radio jet of the 
galaxy 3C303. From this data we derive the magnetoplasma and electrodynamic parameters 
of this 50 kpc long jet. For one component of this jet we obtain for the first time a direct 
determination of a galactic-acale electric current (~ 10^^ A) , and its direction — positive away 
from the AGN. Our analysis strongly supports a model where the jet energy flow is mainly 
electromagnetic. 

Subject headings: galaxies: nuclei — galaxies: magnetic fields — galaxies: jets — quasars: general 



1. Introduction 

A fundamental conundrum of astrophysical jet 
models is how energy is extracted from the accre- 
tion flow close to the black hole event horizon. The 
total energy carried by the jets of active galaxies is 
estimated to be a non-negligible fraction of the su- 
permassive black hole (SMBH) formation energy, 
~ O.lMbhC^ (Kronberg et al. 2001). The jets are 
initially highly relativistic and low-density and for 
this reason are thought to be magnetically domi- 
nated or force-free with a negligible fraction of the 
power in particle kinetic energy. That is, the en- 
ergy outflow from the accretion disk is in the form 
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of a coUimated 'Poynting-flux jet' as proposed by 
Lovelace (1976) and Blandford (1976) and subse- 
quently studied in many papers (Benford 1978; 
Lovelace, Wang, & Sulkanen 1987; Lynden-Bell 
1996; Li, et al. 2001; Lovelace, et al. 2002; 2003; 
Nakamura et al. 2008). The model has a current 
outflow (or inflow) Iz along the spine of the jet of 
cylindrical radius rj initially of the order of the 
Schwarzschild radius of the black hole. The as- 
sociated toroidal magnetic field is responsible for 
collimating the jet. An equal but opposite "re- 
turn current" flows inward (or outward) at much 
larger distances from the jet axis so that the net 
current outflow from the source is zero. Because 
the jet current and its return have opposite signs 
they repel, as a result of their magnetic interaction 
mediated by the toroidal magnetic field B^. This 
repulsion between the jet and its ( more spatially 
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distributed) return eurrent has been demonstrated 
in MHD simulations (Ustyugova et al. 2000, 2006, 
Nakamura et al. 2008). 

Rotation measure gradients observed on parsec 
scales close to the nuclear central black hole have 
provided evidence consistent with an electric cur- 
rent flow (Asada et al. 2002; Gabuzda, Murray, 
& Cronin 2004; Zavala & Taylor 2005). Here, 
on a scale 10* times larger in a "mature" jet, we 
present a measured estimate of the jet's net axial 
current, and its sign. 

2. Brief description of the observations 

3C303 was observed in the VLA's most ex- 
tended "A" configuration in Stokes parameters I, 
Q, U in the 1.4, 4.9 and 15GHz radio bands with 
maximal u — v coverage over an 11.5 hour ob- 
serving period in April 1981. Calibrated images 
have been discussed and analyzed in Kronberg 
(1986) and Lapenta and Kronberg (2005 - here- 
after LK05). In 2010 and 2011, the 1.4 and 4.9 
GHz images were re-edited (flagged), re-calibrated 
and re-imaged using more updated AIPS imag- 
ing procedures, partly clS cl check of the earlier 
1980's imaging procedures. Excellent correspon- 
dence was found, though the more recent pro- 
cedures permitted better imaging of the faintest 
structures. Due to the higher noise levels and 
sparser aperture plane {u — v) coverage of the 15 
GHz data they are not presented here, though they 
are useful for confirming the main image features 
at a resolution of 0.15", and the lack of any 
detectable Faraday rotation between 15 and 4.9 
GHz. 

3. The jet structure and knots of 3C303 

The 3C303 radio source at z = 0.141 (Fig- 
ures 1,2) has a one-sided prominent jet with regu- 
lar repeating knot structures in both VLA radio 
(Kronberg, 1976, 1986; Lapenta & Kronberg 2005 
(LK05), and Chandra X-ray bands (Kataoka et al. 
2003). It resides in a relatively sparsely populated 
intergalactic environment, and at a high galactic 
latitude [{l,b) = -h90°.5, -|-57°.5]. Being away from 
a galaxy cluster, it is well-suited to an analysis of 
Faraday rotation measures (RM) within the jet 
and lobes, since competing Faraday RM from an 
immediate cluster environment and the Galactic 
foreground are small. Given this situation, we ap- 



plied improved, new determinations of the RM's 
of neighboring background radio sources (Kron- 
berg & Newton-McGee 2011). This aUows the 
most accurate available estimate of the "RM zero- 
level", so that our observed RM gradients and 
RM=0 crossings are minimally affected by fore- 
groimd plasma (Section 5). 

The above combination of radio and X-ray im- 
ages plus background source RM's permits the 
jet's physical parameters, including now its cur- 
rent, to be analysed in isolation. 3C303's radiat- 
ing knots are large, a few thousand times larger 
in volume than those in, for example, the well- 
studied nearby M87 radio galaxy. The laterally 
unresolved "spine" of its jet is surrounded by a 
"cocoon" of relativistic, radio and x-ray-visible, 
gas mixed with thermal plasma which has both a 
low plasma - /3 (ratio of plasma pressure to mag- 
netic pressure) and a highly ordered magnetic field 
structure (Fig.l). 

The plasma /? along with the other plasma pa- 
rameter, and minimum energy estimates are de- 
rived and presented in LK05, which can be con- 
sulted for further details. We distinguish the jet 
spine, of cylindrical radius rj normal to the prop- 
agation direction (z) and which is laterally unre- 
solved and appears to carry the jet power, from the 
"knots", which are partially resolved "cocoons" 
around the unresolved "spine" . The cocoons con- 
tain a mix of magnetized thermal and relativistic 
plasma. We find no detectable growth in the knot 
widths away from the AGN, which puts an upper 
limit of 0.7° on the opening angle of the jet. 

The jet's regular structure, combined with the 
appearance of a transverse, magnetically coherent, 
western lobe complex also suggests it as a clear 
laboratory of jet disruption. The disruption point 
is indicated in Figure 1. Model simulations by 
LK05 applied MHD soliton-like solutions of the 
Grad-Shafranov equation to the radio and X-ray 
radio images of the regular knot structure of the 
3C303 jet. Their models successfully computed jet 
stability times and constrained its plasma param- 
eters. The above analysis forms a basis for this 
paper, independent of the MHD model details in 
LK05. It is the visible, synchrotron radiating co- 
coon surrounding the jet that enables us to probe 
the current, most of which is probably confined to 
the unresolved "spine" . It is also apparent that the 
jet continues beyond the disruption point in Fig- 
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lire 1, and terminates in a bow shock- like extrem- 
ity to the west. To the east, a large polarized loop 
is also seen, as well as a possibly related "counter 
hotspot" that aligns with the western jet. We do 
not attempt interepretion of these latter two fea- 
tures in this Letter, and defer them to a following 
paper. 

Below, we describe a rare opportunity to mea- 
sure the Faraday rotation (RM) variation trans- 
verse to, and along, a segment of the jet in knot 
E3 (Figure 1). Combined with the full plasma 
diagnostics, and a determination of the RM zero- 
level from surrounding sources this provides, for 
the first time, a direct estimate of both the mag- 
nitude and sign of a jet's electric current on kilo- 
parsec scales in a low density intergalactic envi- 
ronment. 

The partial resolution of the knots in the VLA 
images at 1.4 and 4.9 GHz permits estimates of the 
relativistic gas density and approximate magnetic 
field energy within each kpc— scale synchrotron 
radiating knot (cocoon). 

The first stage in measuring the jet current is to 
calculate the approximate magnetic field strengths 
in the knots using the analysis described in LK05. 
Note that the knots' cquipartition magnetic field 
strength estimates depend only weakly on the vol- 
ume filling (/), as , and on the upper cutoff 
frequency of the photon spectrum of the jet's emis- 
sion. The latter extends to 10^^ Hz (Kataoka et 
al, 2003), and its reduction by, say xlOO, would 
reduce the (rcdshift-corrected) total energy con- 
tent in the knots by only 30%. At lower particle 
energies the uncertain extrapolation of the knots' 
relativistic particle energy spectrum, correspond- 
ing to (unobserved) very low radio frequencies, 
could in principle raise our estimate of \B\ by a 
modest factor less than unity, however no anoma- 
lously steep spectrum between 1.4 and 4.9 GHz is 
seen in the region of the 3C303 jet. 

We estimate a field strength in 3C303's jet 
knots, where Bk^ot ~ 0.5 milligauss (LK05), mea- 
sured with a 0.35" beam and allowing for some 
out-of-plane component due to the expected he- 
lical form of the jet's magnetic field. 

At the 0.35" resolution of the 4.9GHz VLA im- 
age (Figure 2) the projected line-of-sight magnetic 
field orientation within each knot is parallel to the 
jet axis. The 4.9GHz percentage polarizations of 



knots El, E2, and E3 are - 27%, 21%, and 26%, 
respectively, and the knots are elongated in the 
projected jet direction by ~ 2.2". 

4. The Faraday rotation and thermal elec- 
tron density 

Figure 3 shows the 2-D RM variation within the 
brighter parts of the entire 3C303 system derived 
from the 1.4 and 4.9 GHz images at a common 1.5" 
resolution. Values range between -|-5 and -|-25 rad 
over a maximum extent of ^ 50", including 
3C303's extended western lobe. No higher RM's 
are seen anywhere within the source boundaries. 
One test for a 180° ambiguity in xW if' to sep- 
arately produce Q and U images at the higher 
frequency of 5GHz, where we can test for "ex- 
changes" between Q and U along the jet. None 
were found, meaning that no substantially higher 
RM's occured in or around the jet. The ambiguity- 
free RM of 3C303's integrated emission (from po- 
larization data at 6 frequencies) is -|- 18 ± 2 rad 
m~^ (Simard-Normandin et al. 1981). This is 
nicely concordant with the RM range in Figure 
3, after the -18 rad m~2 correction above for the 
independent foreground RM estimate in §5 below. 

The jet knot E3 is resolved just enough both 
across and along the jet, to show an RM gradient 
which is independently transverse to the jet axis. 
The varying RM in knot E3 also independently 
confirms the presence of some thermal plasma in 
the cocoon. The RM gradient, VRM around 
knot E3 is ~ 10 radians m~^kpc~^. Because the 
RM in Fig. 3 is resolved in 2D over 10 - 12 inde- 
pendent sampling points in a projected y — z plane 
around knot E3, we are able to specify both the 
magnitude and direction of WRM within the knot 
E3 cocoon. 

Having measured a differential RM over the 
knot, its dimension, and the approximate mag- 
netic field strength, we can estimate the non- 
relativistic plasma density within the knot via the 
following equation. 




where x is the line-of-sight distance through the 
knot. Inserting B|| and x gives Ug ~ 1.4 x 
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10 ^cm similar to the value obtained in LK05. 

An independent estimate of rig can be made 
from the observed polarization degree in the knots. 
On the reasonable assumption that the very regu- 
lar field geometry in knot E3 results from a helical 
geometry, knot-internal Faraday dispersion is lim- 
ited by the high observed degree of polarization 
~ 25% at 5 GHz. This test implies a back-to- 
front, knot-internal differential Faraday rotation 
at 5GHz of <^ 30°, consistent with the knot's RM 
observed above (LK05). It is consistent with the 
above rig estimate and gives us confidence in the 
local Ue scaling that is required to estimate the jet 
current - which we discuss next. 

This low value of rig outside the jet spine points 
to a magnetically confined (or dominated) jet 
(LK05): Confinement of the jet by hot gas at X- 
ray temperatures is not possible because pgas = 
rie kTx < bVStt. a common feature of the mag- 
netic jets is that the internal pressure of the mag- 
netic field B^jjq^/Stt (and relativistic particles) is 
confined by the pinch force of the external toroidal 
magnetic field The models give -Bknot ^ -By 
(LK05). 

5. Estimate of the electric current in the 
3C303 jet 

An important final step in relating (VRM) to 
an axial current is to calibrate the RM zero level. 
This was done by averaging discrete source RM's 
along neighboring lines of sight to 3C303 from a re- 
cent accurate RM dataset (Simard-Normandin et 
al. 1981, Kronberg & Newton-McGee 2011). The 
resulting RM zero-level correction, —18 ±4 rad 
m~^, has been applied to the RM scale of Figure 
3. With this zero- level correction, we find that 
the RM in knot E3 changes sign on the jet axis, 
within the limit of our uncertainties. Furthermore, 
its gradient (VRM = n^B^) is found to be perpen- 
dicular to the (independently measured) axis of 
knot E3 and the jet axis. This corresponds to the 
expected RM behavior around an electric current 
flowing along a jet axis. 

The reversal of the RM sign on the jet axis 
matches the signature of a current-carrying cylin- 
der, where we are measuring the azimuthal com- 
ponent of a helical magnetic field at distances 
« y on opposite sides of the synchrotron radiat- 
ing cocoon surrounding the unresolved spine of the 



jet. In terms the measurables, the inferred cur- 
rent, Iz oc {y/ne)WRM is the total current within 
the cocoon's radius y. Over this distance we mea- 
sure VRM of ^ -I- 10 rad m'^/kpc. Thus, within 
knot E3 

where S,^ is the toroidal component of magnetic 
field at a distance y in the sky plane from the jet 
(z) axis. The positive sign of the gradient V(RM) 
indicates that a current of ~3.85 x 10^* A, for 
Sequi=0.5 mG, is directed away from the galaxy 
nucleus. That is, within 0.5 kpc of the jet axis 
there is a net flow of negative charge toward the 
AGN core within knot E3. 

Our result for the jet current is subject to differ- 
ent corrections and adjustments: For example, (1) 
if the jet axis is significantly out of the plane of the 
sky then the observed RM gradient can arise from 
a purely toroidal field, but a helical field would 
in this case shift the RM= line away from the 
jet's axis, contrary to the present observations. (2) 
In another example if the magnetic field of the 
jet is not symmetric about its axis the field esti- 
mates will be modified. Such possibilities, espe- 
cially (1), underline the importance of establish- 
ing from observations, an optimally foreground- 
free RM {y,z) over the source image, as wc have 
done here. Considering these uncertainties and 
the limited angular resolution, our best estimate 
of the axial jet current at knot E3 is iQi^-^^O-^ 
amperes. This is the first direct estimate of a cur- 
rent in a kpc-scale extragalactic jet. 

6. Discussion and Implications 

A magnetically dominated jet can be modeled 
as a transmission line running in the direction 
with the electric potential drop across it AF. This 
is equal to the potential drop across the BH accre- 
tion disk from its inner- to outer radius (Lovelace 
1976; Lovelace and Ruchti 1983). A typical value 
is AV ~ lO2OV(Sbh/lO'*G)(Mbh/lO«M0), where 
is the poloidal magnetic field strength near the 
black hole and Mbh is its mass (Lovelace 1976). 
The current flow in this "transmission line" is 
Iz = /SVjZ, where the impedance of a relativistic 
jet isZ ^ c-i(cgs) = (4^)-i(^o/eo)'/'(MKS) = 
30fl (Lovelace 1976). More generally, for a jet with 
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bulk axial velocity u^, Z — {l/c){tiz/c), (Lovelace 
& Romanova 2003). The electromagnetic power 
transported by the jet is Ljot — I^Z . 

We can now apply these ideas to the 3C303 sys- 
tem. The total photon luminosity Lrad integrated 
from lO^Hz to lO^^Hz is about 3.7 x 10^^ erg s'^ 
or 3.7x lO^^W (LK05). The radiated power is ex- 
pected to be significantly less than the jet power, 
because the jet power goes into a combination of 
pdV work on the ambient medium, and to energiz- 
ing electrons and ions around the jet and in the 
outer radiating lobes. Thus, we let L^ad = e^jet, 
where £ < 1 is an cifficicmcy factor. Substituting 
the observed values of Lrad and 7^ we find 



Tvrad 



(3) 
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-3 



' £rad \ / 3.85X IQiSA y 



0.1c" 



where we have normalized the jet bulk velocity to 
0.1c. Note that e is safely less than unity, and the 
model of Poynting flux transport of jet energy in 
3C 303 's black hole/accretion disk system appears 
to be self-consistent. 

The system evidently needs to incorporate a 
"transducer" that converts the Poynting energy 
flux into high energy particles which then pro- 
duce synchrotron radiation. These issues, and 
the complex lobe structure will be discussed in a 
subsequent paper. 
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Fig. 1. — 1.4 GHz, 1.5" resolution VLA image 
of the entire 3C303 system showing an appar- 
ent counter-jet knot, the highly polarized west- 
ern lobe, and a continuation of the jet to its final 
western stopping point. Linear polarization lines 
are scaled such that a 1" equivalent length = 2.5 
mjy/beam. I -contours are shown at -0.5, 0.5, 1, 
1.5, 2, 3, 6, 12, 24, 50, 100, and 150 mJy/beam. 
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3C303 4866 MHz 
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Fig. 2.— A 4.9GHz VLA image of the 3C303 jet 
at 0.35" resolution showing the 3 prominent, elon- 
gated and equally spaced knots El, E2, and E3 to 
the right of the stronger, variable milliarcsec-size 
galaxy nucleus source. 
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Fig. 3. — Faraday rotation image of the 3C303 
radio source at a resolution of 1.5". The RM 
zero level has been corrected by 18 ±4 rad ni~^ 
(see text). 
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